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We present an analysis of the dynamics of the equifacial pentahedron on the Kapovich-Millson 
phase space under a volume preserving Hamiltonian. The classical dynamics of polyhedra under 
such a Hamiltonian may arise from the classical limit of the node volume operators in loop quantum 
gravity. The pentahedron is the simplest nontrivial polyhedron for which the dynamics may be 
chaotic. We consider the distribution of polyhedral configurations throughout the space and find 
indications that the borders between certain configurations act as separatrices. We examine the 
local stability of trajectories within this phase space and find that locally unstable regions domi- 
nate although extended stable regions are present. Canonical and microcanonical estimates of the 
Kolmogorov-Sinai entropy suggest that the pentahedron is a strongly chaotic system. The pres- 
ence of chaos is further suggested by calculations of intermediate time Lyapunov exponents which 
saturate to non zero values. 



I. INTRODUCTION 

It has long been known that black holes act as thermodynamic systems whose entropy is proportional to the area 
of their horizon [I] and a temperature that is inversely proportional to their mass [2]. This has raised the question 
how quickly and how perfectly a black hole thermalizes and thereby effectively destroys any injected information. 
Sekino and Susskind have argued that black holes are "fast scramblers" [311], i. e. systems that render information 
practically unretrievable at the maximal possible rate. 

Scrambling information among different degrees of freedom is one particular form of deterministic chaos. A well 
known realization of scrambling is provided by the Baker's map. Other examples of dynamical chaos generated by 
the nonlinearity of Einstein's field equations have been studied in the context of various analytical solutions of general 
relativity (see [3 E] for a review). The results of these investigations suggest that dynamical chaos, and thus the 
tendency to lose initially known information effectively irretrievably, is a generic property of classical gravitation. 

These considerations apply to macroscopic gravitational fields, which can be considered as thermodynamic systems 
with a macroscopic number of degrees of freedom. But what about microscopic black holes with masses near the 
Planck mass, which possess only a small number of degrees of freedom? Is there a smallest black hole that can be 
considered as a thermal system? Which mechanism drives the apparent thermal equilibration of black holes at the 
microscopic level? The pursuit of these questions requires a quantum theory of gravity. While the true quantum 
mechanical foundation of gravity is still unknown, there are at least two widely explored candidates for such a theory: 
superstring theory [3 [8] and Loop Quantum Gravity [5Ull] . Here we consider the problem of the microscopic origin 
of the thermal properties of space-time in the framework of Loop Quantum Gravity (LQG) . 

LQG is an attempt to reconcile general relativity and quantum field theory [TfJl [T^HT^] . the structure of space-time 
emerges naturally from the dynamics of a graph of SU (2) spins. Naively the nodes of this graph can be thought of as 
representing granules of space-time, the spins connecting these nodes can be thought of as the faces of these granules. 
The volume of these granules, along with the areas of the connected faces are naturally quantized |15j . 

There has recently been a focus on finding a semi-classical description of the spectrum of the volume operator 
at one of these nodes. There have been several reasonable candidates for the quantum volume operator (cites), a 
semi-classical limit may serve to pick out a particular one of these forms. The volume preserving deformation of 
polyhedra has recently emerged as a candidate for this semi-classical limit [151 HZ] ■ I n this scheme the black hole 
thermodynamics can be derived [TB] in the limit of a large number N of polyhedral faces. Here the deformation 
dynamics of the polyhedron is a secondary contribution after the configuration entropy of the polyhedron, which can 
be readily developed from the statistical mechanics of polymers. 

The dynamics of the elementary polyhedron, the tetrahedron, can be exactly solved and semi-classically quantized 
through the Bohr-Sommerfeld procedure |19l 120] . The volume spectrum arising from quantizing this classical system 
has shown remarkable agreement with fully quantum calculations. If the tetrahedron is the "hydrogen atom" of space, 
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the next complex polyhedron, the pentahedron (N = 5), can be considered as the analogue of the helium atom. Just 
as the full range of atomic physics phenomena occurs first in the helium atoms with its two electrons, the pentahedron 
represents the first space configuration that puts the ideas discussed above to a non-trivial test. 

The dynamical system corresponding to the isochoric pentahedron with fixed face areas has a four-dimensional phase 
space compared with two dimensional phase space of the tetrahedron. The integrability of a Hamiltonian system in 
a phase space with more than two dimensions is not assured. Non-integrable Hamiltonian systems exhibit a wide 
variety of interesting behaviors, including Hamiltonian chaos. The incompressibility of Hamiltonian flows leads to a 
very specific form of chaos j2 1-25] . A well known example is the three-body Kepler problem, which is non-integrable, 
in contrast to the two-body problem [21] and exhibits chaos at the classical level. In the analogous quantum system, 
the helium atom, the chaotic dynamics exhibits itself through the mixing of bound and continuum states. 

There are two distinct classes of polyhedra with five faces, the triangular prism and a pyramid with a quadrilateral 
base. The latter forms a measure zero subset of allowed configurations as its construction requires reducing one of 
the edges of the triangular prism to zero length. This process imposes an algebraic constraint between configuration 
variables making this a co-dimension 1 configuration. 

In this article we first briefly review the symplectic Kapovich-Millson phase space of polyhedral configurations in 
Section [il] In Section III we outline a method by which it is possible to uniquely construct a triangular prism (or 



quadrilateral pyramid) for each point in the four-dimensional phase space and review a method due to Laserre (25j for 
computing the volume of any polyhedron from its face areas and their normals. We introduce a convenient labeling 
convention of the faces of the pentahedron in Section |IV| With these tools in hand we are then able to compute the 
trajectory of a pentahedron starting from a given point in phase space with the volume of the pentahedron serving as 
its Hamiltonian. This volume preserving (isochoric) evolution keeps the areas of the faces of the pentahedron fixed 
and only allows for changes in their shape and orientation. In this way the edge lengths and vertex positions of the 
pentahedron are free to vary, in effect the pentahedron is smoothly deformed, keeping its volume constant. In Section 
[V] we examine the local dynamical stability of configurations of the pentahedron throughout the phase space and 
map the distribution of the local Lyapunov exponents as function of the volume of the pentahedron. To explore the 
stability further we compute the time evolution of intermediate time Lyapunov exponents which suggest the presence 
of chaos at small volumes. We conclude our article with a summary of our results and a brief analysis of their possible 



implications in Section VI 



II. POLYHEDRA AND PHASE SPACE 

A convex polyhedron is a collection of faces bounded with any number of vertices. A theorem by Minkowski |26] 
states that the areas A\ and normals ni of each face are sufficient to uniquely characterize a polyhedron. If we define 
A\ = Aifti then the polyhedral closure relationship 

i 

is a sufficient condition on Ai to uniquely define a polyhedron with N faces. The space of shapes of polyhedra Vn 
with N faces was further investigated by Kapovich and Millson [27] • The shape space is defined as the space of all 
polyhedra modulo to their orientation in three-dimensional space: 

V N = ^A l \Y^A l =Q,\A l \=A^/SO{?,). (2) 

The shape space of convex polyhedra with N faces is thus 2(A^ — 3)-dimensional; in particular, the shape space of the 
tetrahedron (N = 4) is two-dimensional and that of the pentahedron (N = 5) is four-dimensional. As Kapovich and 
Millson showed, this space admits a symplectic structure, which can be defined by introducing a Poisson bracket for 
any two functions f{A{),g(A{) as 

Canonical variables with respect to this Poisson bracket are defined on Vn as follows. One first defines the vector 
sum pk = Yli=i Ai of the first k + 1 oriented faces. The ordering/labeling of these normals is not important for the 
physical shape of the polyhedron. Minkowski's theorem ensures that any set of vectors which obey (JlJ will produce 
a unique polyhedron. 
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The canonical momenta in the KM space are defined as pk — \pk\ and the conjugate positions are given by the 
angle between pk X Ak+i and pk x A^+2- Using ^ one can verify that {qk,p'k} = ^fcfc'- We refer the reader to [T51I27] 
for detailed explanations and proofs. 

These quantities may be visualized by representing the polyhedron as a polygon with edges given by the vectors 
Vi = Ya=i M- Note that this generally yields a non-planar polygon. Now consider systematically triangulating this 
polygon starting from one vertex, which is chosen as the origin. The inserted edges are the conjugate momenta p and 
the angles between each of these edges are the conjugate positions. An illustration of the pentagon associated with a 
pentahedron in shown in Fig.JT] 

The momentum vectors pk define axes about which subsets of the normals rotate as the polyhedron is deformed; 
the angles qk measure the rotation about these axes. The collection of coordinates (pk,qk) thus completely describes 
the "bending flow" of the polyhedron. 

The momentum vectors pi,p2, ■ ■ ■ ,pn partition the polygon into a set of coupled subsystems. In Fig. [2] we show 
several pentagons illustrating the variation of the angles q\ (left) and q2 (right), while holding the remaining phase 
space variables fixed. 




FIG. 1: An example configuration of the system in the polygon representation (left), the phase space coordinates plotted here 
are z = {0.3,0.4,0.9,0.91}. The normal vectors are plotted as the red solid arrows and the momentum vectors are plotted as 
the dashed blue arrows. We also show a rendering of the associated polyhedron (right). All polyhedral faces have area fixed to 
one, so all polygonal edges have unit length . 

It is not trivial to calculate the volume of a general polyhedron directly from the phase space coordinates. Haggard 
has computed an analytic formulation for the |2S]. Alternatively it is straightforward, if cumbersome, to compute the 
volume of a polyhedron from its normals and face areas. We seek a mapping F : z — > n{z) between a point z in the 
phase space and the normals n(z) associated with that point. 

We now specialize our discussion to the pentahedron, although it is easily generalized to higher polyhedra. One of 
the normals can be completely defined in terms of the four others by the closure of the pentahedron. Since we require 
that the normal vectors be normalized, so that the areas of the faces are fixed, this leaves two degrees of freedom per 
vector. After enforcing closure we thus need seven equations to completely specify the pentahedron in space. The 
definitions of the canonical variables provide four equations, leaving three free quantities. We fix these three final 
components by picking an orientation of the coordinate system. This is sufficient to specify F : z —¥ n{z). 

It is convenient to orient the polyhedron in such a way that one of the faces lies in the x—y plane and its normal is 
oriented in the negative z direction (0, 0, —1). To fix the final coordinate direction we require that one of the edges of 
this face is parallel to the x direction. This fixing the normal of the adjacent face to have components (0, n y ,n z ). We 
are free to label our normals in any order. Let n 5 be the normal fixed by closure, ri\ be the normal in the x—y plane 
and 712 the normal to the face which has the edge parallel to x. 

It is useful to explicitly derive the inverse mapping F^ 1 : n(z) — > z, i. e. the map between a pentahedron represented 
by a set of normals and face areas and the set of canonical variables z = {q%, q2,p\,p2\- The canonical variables can 



FIG. 2: A superposition of configurations showing the rotation of qi (left) and 52 (right) through n/2 and 7r radians. The 
remaining coordinates are held fixed. The normal vectors are plotted as the red solid arrows and the momentum vectors are 
plotted as the dashed blue arrows. All faces have area fixed to one. 



be obtained from the normals and areas as: 
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It is helpful to introduce the following vectors 



pi = \Ai +A 2 \ = \-A 5 -A A - A 3 \, 
P2 = |Ai + A 2 + A 3 1 = I - A 5 ~ A 4 \, 
q 1 = angle {pi x A 2 ,px x 1 3 \ , 

angle W 2 X A 3l p 2 x AA ■ (4) 



Po = Ax, 

pi = Ai + A 2 = p + A 2 , 

p 2 = A x + A 2 + l a = pi + 1 3 , 

p 3 = Ax + A 2 +1 3 +1 4 =p 2 + 1 4 , (5) 



the canonical angles can be written as 



qx = angle {p Xpx,px X p 2 } , 

q 2 = angle {pi xp 2 ,p 2 x p 3 } . (6) 

To ensure that the angles qx , q 2 are uniquely defined over the phase space we take care to compute the signed angles 
between the vectors po x pi and pi x p 2 . We define the signed angle between two vectors a, b with a reference vector 
r in ([7]). In practice the reference vector is the appropriate momentum vector p^- This is sufficient to define F~ x . 



S 

c 

a 



I a x b I 

a ■ b 
r-^axbj, 

I arctan(C, S) if sgna > 0, ,^ 
I 27r — arctan(C, S) if sgn a < 
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Given the choice of orientation of our coordinate system along with face areas A^,,,^, the polyhedral configuration is 
uniquely determined by the following five vectors 

n x = (0,0,-1), 
n 2 = (Q,n 2 y,n 2z ), 
n.3 = (n 3x ,n 3y ,n 3z ), 
n± = (ri4x, n Ay , n 4z ), 

n 5 = --j- (Aifii + A 2 n 2 + A 3 n 3 + A 4 n 4 ) . (8) 

^5 

We use normalization to fix the magnitudes of n 2y ,n 3x n 4xi leaving 5 unknown components. The following set of 
equations arise from inserting Q into the definitions of \p\ \,\p 2 \ and \p 2 \ — \p\ + Ax \ , 

p? = A\ + A\ + 2A l A 2 n 1 ■ n 2 , (9) 
pj = Al + Al + 2A 4 A 5 n 4 -n 5 , (10) 
pi = p\ + A\ + 2A 2 A 3 n 2 ■ n 3 + 2A 1 A 3 n 1 ■ n 3 , (11) 

To close the system we use the definitions of the canonical angles ffify. At this point the system can be inverted 
numerically to give solutions for the components of the normals in terms of the phase space variables. Particular care 
needs to be taken to ensure that unique configurations of normals are obtained for each set of phase space variables 
z. Further care needs to be taken in picking a sign for the components whose magnitudes are fixed by normalization. 
The signs of these components vary across the space in a non trivial fashion. We introduced a set of heuristics to 
suggest which set of signs might be appropriate, the inverse map F -1 is vital for testing and rejecting trial solutions. 

It is helpful to note that the set {Ai,A 2 ,A 3 ,—p 2 } can be taken as the faces of a tetrahedron with canonical 
momentum p\ and conjugate angle q\. The dynamics of this tetrahedron are independent of the second bending angle 
q 2 . The set of equations to be inverted can be factored into a set of 3 equations for n 2zi n 3y ,n 3z in Pi,p 2 ,qi and a 
set of 2 equations for n 4y , n 4z which are dependent upon the full set of canonical variables. The former can be solved 
analytically. Fixing all face areas to A and restricting q\ to the first quadrant gives 



= ( 12 ) 
= -(4A 2 -pj)(Ai +pj -pi) ± |cos(g 1 )| v /(-4^ + p\) (g + (g _ p 2 )2 _ 2 ^(pf 

4A 2 Pl ,/AA 2 - p\ 



^ ri v n 

(A 2 +pl -p\) ± |cos( gi )|^(-4A2 +p\)(A± + {pl-plY - 2A>{p\+pl))^ . (14) 

The full solution with different face areas and with q± taking values in all quadrants is obtained using a computer alge- 
bra system. A Newton based root finder [21] is then used to find simultaneous solutions to q 2 = angle {p 4 x p 2 ,p 2 x p 3 } 



U3z ~ AA* Pl 



and ( 10 1 



The numerical integration of the dynamical system, which is based around gradients of the Hamiltonian, requires 
considerable numerical accuracy in F. The choice of initial guess for the root search when solving for n 4yy n 4z is 
vital for obtaining a solution. To improve this process an empirical PDF with density inversely proportional to the 
euclidean norm of the residual vector from the root finding routine was generated by evaluating trial solutions over a 
grid in the space. This PDF was simultaneously refined while being used to generate initial guesses for the roots. 

The construction method outlined above can be generalized to develop a chain of polynomials for N > 5 allowing 
at least numerical evaluation of F. 



A. The shape of the phase space 



The geometric structure of the polyhedron itself, particularly the fixed face areas, induces certain restrictions upon 
the phase space. The position space is 2tt periodic by construction. The momentum space is restricted by the areas 
of the faces, from the triangle inequality 

Pi = |i , i+i 2 |<|i'i| + |i!a|, 
P2 = \pi+A 3 \ < \pi\ + \A 3 \- 
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Heron's formula for the area of a triangle can be used to simplify the above inequalities, 



2t(a,6,c) = -p/{a + b + c)(a + b - c)(a - b + c)(b + c - a) 



(15) 



where a,b,c are the edges of the triangle. Consider the triangles Ai = {Ai, A2,pi}, A2 = {pi, A^,^} and A3 = 
{p2,A4, A 5 }. For the system to be in a reasonable configuration we require that the area of each of these triangles be 
non zero, i.e that the terms under the radical in ( 15 ) be positive. This restriction on Ai and A3 implies 



Applying this to A 2 gives 



\A 1 -A 2 \ < Pl < \A 1 +A 2 \ 
\A 4 -A 5 \ <p 2 < \A 4 + A 5 \ 



M + Pl > P2, 
A3 +P2> Pi, 
P1+P2 > A3, 
+ Pi +P2 > 0. 



(16) 



(17) 



The analysis in the remainder of this article is focused on the case where all the areas are fixed equal to one, the 
allowed momentum space for this case is shown in Fig. [3] 
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FIG. 3: The shaded area shows the allowed region of momentum space, satisfying ( 16 1 and (171, for a system with Ai — A2 
A 3 = A 4 = A 5 = 1 



III. HAMILTONIAN & POLYHEDRAL RECONSTRUCTION 



Following from [T8H2T)] we use the volume of the pentahedron at a given point in the phase space as the Hamiltonian. 
This ensures that trajectories generated by Hamilton's equations will deform the pentahedron while maintaining a 
constant volume. Unlike the case of the tetrahedron, it is not trivial to compute the volume of the general pentahedron 
in terms of the canonical coordinates. Although elegant direct geometric expressions for the volume of any polyhedron 
exist |28j . these do not give particularly tractable expressions in terms of the phase space variables. 

Consider a vector field F(x) — using the divergence theorem we can find the volume of a polyhedron 

V = [ V • Fdfl = I F ■ nds = V -x*i • Hi (18) 
Jq Js 3 

where f2 is the interior of the polyhedron, is a point on the i'th face and rij is the normal to that face. 



We can compute the volume of a polyhedron specified as a set of normals and areas using ( 18 ) once we know the 
location of a point upon each face. To obtain this we need the edges, or equivalently vertices, of the polyhedron. 
These can be found by reconstructing the polyhedron from its normals and areas, a process originally due to Laserre 
[25j and more recently outlined in |18j . This Laserre reconstruction procedure works for any N and requires only the 
minimization of a quadratic function, in some cases algebraic results can be directly obtained. 



FIG. 4: A set of renderings of the reconstructed pentahedron corresponding to the point z\ = {0.3,0.4,0.9,0.91} in the KM 
space, all face areas set are to 1. Where visible internal edges have been drawn with dashed lines, external edges with solid 
lines. 



The reconstruction process works by pushing a set of N infinite planes, some set of distances h = {hi, . . . , hpj} from 
the origin along the normals of the system. The areas of the polyhedral faces formed by the union of the half-spaces 
defined by these planes can be computed as a function of the heights A(h). A numerical minimization routine is 
then used to extract the set of heights which minimize the form 

\A(h)-A\ 2 

where A is the set of desired areas. (In passing we note that by carrying out this procedure after flipping the sign of 
the normals we can obtain the chiral dual of a given polyhedron). The set of vectors {/iini, . . . , /i/v^jv} can then be 



used in ( 18 ) to compute the volume of the system. 

As an example we reconstruct the pentahedron corresponding to z\ = {0.3,0.4,0.9,0.91} shown in Fig. [I] with all 
the face areas fixed to 1. Applying our mapping to compute the normals F(z\) gives: 

ni = {0,0,-1}, 

n 2 = {0,0.803,0.595}, 

n 3 = {0.249,-0.114,0.961}, 

n A = {0.301,-0.921,0.244}, 

n 5 = {-0.550, 0.232, -0.801}. (19) 

After applying the reconstruction routine we obtain the following set of heights 

h = {0.237, 0.239, 0.236, 0.239, 0.237} (20) 

with a corresponding volume of 0.396621. The heights found are of a similar magnitude because the face areas are 
all fixed to the same value, this precludes very anisotropic configurations. A set of renderings of the polyhedron is 
shown in Fig. [3] 

In the remainder of this article we will consider the case where all faces of the polyhedron have the unit area. We 
can estimate some limiting values of the volume in this case. A regular equilateral triangular prism should be close 
to the maximum volume. If the triangular edges have length I and the vertical quadrilateral edges have length h, 



the volume is Vprism = A£, where A denotes the common area of all faces. Using ( 15 ) we obtain the triangular edge 
length I = 1^2, and the quadrilateral edge length h = A/£ so V plism = ^A 3 / 2 sa 0.658A 3 / 2 . The volume of a 
regular square pyramid can also be obtained from V pyr = \Ah where A is the area of the base and h is the height of 
the pyramid. Again using (15) we obtain h = \/l5A/2 and so U pyr — ^^-A 3 ^ 2 w 0.645A 3 / 2 . It is interesting to note 



that the regular square pyramid has a slightly smaller volume than the regular triangular prism. 

The volume is bounded from below by zero, in the limit of collinear normals the volume is certain to vanish. A 
general feature of this Minkowski polyhedral reconstruction process is that if the area of one face is smoothly shrunk 
to zero the system will smoothly deform into a polyhedron of a lower order |18j . 
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FIG. 5: A section in the qi,<Z2 plane through the Hamiltonian evaluated at p 1 = p2 = 0.94, all face areas are fixed to f. The 
contours are isochors, the color scheme is brighter at larger volumes. 



We can estimate the behavior of the volume as it approaches zero by considering the extreme cases of ( |I7[ ). As we 
shrink the second momentum p 2 to zero then p% — > A. In this limit the area of the triangle A3 = {A4, A$,p2\ goes to 
zero, we suppose that the volume of the system is then entirely dominated by the set of vectors {Ai, A 2 , A 3 ,p 2 }, the 
volume of the associated tetrahedron []j|J [50] is 



a 8%(A 1 ,A 2 ,p 1 )%(pi,p2,A 3 ) . , . 
V tet = g sm( 9l ), 



(21) 



for a fixed value of p\ on the boundary of the allowed space the volume goes smoothly to zero as p 2 is taken to its 
extreme value. Taking a unit area system for simplicity with q\ = n/2,p± = 1,P2 = 1 — % and expanding (21) for 
small x we obtain 



Vi 



tet\A=l 
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The same construction can be carried out mutatis mutandis to obtain a limiting a form of the volume as the first 
momentum is shrunk to zero. The derivatives of the volume along the limiting direction are well behaved in this limit. 

A section of the polyhedral volume through the q\ , q 2 plane is shown in Fig.[5j the spots are regions where the 
numerical reconstruction routine failed. The contours show lines of constant volume (isochors), cooler colored regions 
correspond to smaller volumes. The volume is visibly 27r periodic in both position variables. The structure of the 
Hamiltonian in this projection suggests that of a complicated coupled harmonic oscillator in a periodic space. In Fig. [6] 
we show a projection of the Hamiltonian in the pi,p 2 plane at fixed qi,q 2 . The lowest volume regions in this plot are 
found where the momenta take their minimum allowed values. The dependence of the volume upon the momenta is 
not trivial, the Hamiltonian appears to be non separable in its current form. 

In Fig. 14 we show a set of projections in the pi,P2 plane which together span the angle space from ir to 2ir. The 



behavior in the other quadrants of the q\ , q 2 range is very similar to the one shown here, we present this reduced set 




FIG. 6: A section in the pi,P2 plane through the Hamiltonian evaluated at q\ = 0.3, 52 = 0.4, all face areas are fixed to 1. The 
contours are isochors , the color scheme is brighter at larger volumes. 



9 



of plots to give the reader some intuition of the variation in the volume over the whole space. The scarring visible in 
some of the subplots is a result of the failure of the root finding heuristic along boundaries where m x flips sign. The 
volume is symmetric about the line q\ = Q2 and the maximum volume is modulated sinusoidally as the angles pass 
through a period. The maximum volume is seen to be on the p\ — P2 line for q± = qi ~ Btt/2, in the other angular 
quadrants the maximal volume is obtained at multiples of tt/2. 

The plots around maximum volume are symmetric about the line p\ = P2, here the p\,P2 projection is well described 
by the sum of the areas of the triangles A 1; A 2 , A 3 as given by (Tl5| 



. vSinfoil , , . sinfoi ) sinfgo) , . , sulfa?) 

H( qi ,q 2 , Pl ,p 2 ) ocf&fa, A lt A 2 )^^ l +%(pi,P2,A 3 ) p> f> + Si (p 2 , A 4 , A 5 )—^, 22 



this is a good description across the whole range of momenta and for small angular deviations around 7r/2. 



This choice of the form of ( 22 ) is motivated by the expression for the squared volume of a tetrahedron obtained in 
[T5] , The radicals in the denominator are needed to fix the dimensions. The canonical momenta have dimensions of 
area so the Hamiltonian must have dimensions of A 3 / 2 . The full Hamiltonian is deformed in a more complex way 



when moving away from multiples of ir/2 in the angle space. The approximate Hamiltonian (22 1 may give insight 
into the dynamics of the system in the large volume limit, we shall present an analysis of the trajectories of this and 
the full Hamiltonian in a future work. 

An extended set of projections in the qi,q 2 plane is shown in Fig. |15| the set of plots shown spans the full range 
of both momentum variables. The periodic structure shown in Fig. [5] exists throughout the whole space, however the 
periodicity is carried by volume maxima at small momenta (bottom left) and volume minima at large momenta (top 
right). 

The volume of the triangular prism has a very complex structure over the allowed regions of the KM phase space. 
The volume is observed to be smooth and differentiable across the space. Numerical examination of the volume in the 
limit that the system approaches the momentum boundaries of the phase space indicate a smooth approach to zero. 



IV. CONFIGURATIONS 

Due to our particular choice of the orientation of the R 3 coordinate axes we are able to uniquely distinguish the 
normals in our implementation of the triangular prism. Our numerical map from a phase space point z to the normals 
n(z) makes identification of each of the remaining normals clear and gives them a natural ordering. 

We adopt a similar scheme to Haggard [28] and label the possible configurations of the pentahedron in terms of 
which normals correspond to square or triangular faces. Note that a full reconstruction of the pentahedron is required 
to determine the number of edges (or order) associated with each face and therefore with each normal. Fixing all the 
faces to have the same area we are no longer able to uniquely identify a particular square or triangular face. However 
we are still able to distinguish the order of the faces. 

We label the distinct configurations by concatenating the indices of the triangular faces. A configuration where the 
third and fourth faces are triangular would be labeled 34 or equivalently 43. We label the configurations in lexical 
order to remove this ambiguity. The set of available configurations C is 

C = {12, 13, 14, 15, 23, 24, 25, 34, 35, 45}. 

The triangular prism can be smoothly deformed from one configuration to several others. These deformations 
involve a two step process. First the smooth shrinking of one the edges shared by two square faces to zero, creating a 
quadrilateral pyramid in the limit that the edge vanishes. The vertex at the apex of this new pyramid is split into an 
edge which now borders two newly square faces. An example of this process in shown in Fig. [7] here a 13 triangular 
prism configuration is deformed into a 24. 




FIG. 7: Schlegel diagrams showing the transformation of a 13 configuration to a 24 configuration. The 5 face, not shown in 
these figures, forms the back face of the graphs. 
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The undirected graph in Fig. [8] shows the set of allowed deformations. For example the transformation 13 — > 24 
shown in Fig. [7] would be represented as moving from the blue vertex to the dark green vertex. The colors of the 
vertices in the graph serves as a key to the colors shown in Fig. [16] and Fig. [17] These figures show the configuration 
of the system across the phase space in the pi,P2 and qi,q2 planes. The graph can be used to understand which 
boundaries in these figures can be smoothly crossed and which are disallowed. 




24 14 



FIG. 8: Allowed transitions between configurations. The vertex colors correspond to the colors plotted in the configuration 
scan figures Fig.[l6] Fig.[T7] 

The set of momentum-space projections shown in Fig. [16] show a complex but smooth distribution of the interfaces 
between configurations through the space. Many of the boundaries represent allowed transitions. Several disallowed 
transitions are also observed. The 34 configuration (red) is spatially adjacent to the 45 configurations (light green) in 
the first column of figures, in the first row the 23 configuration (yellow) is adjacent to the 12 (gray). These unphysical 
transitions may serve as a kind of separatrix in the momentum space. The position-space projections shown in Fig. |17| 
exhibit a greater complexity in the boundaries between the configurations. However disallowed transitions are not 
obvious in this projection. 

V. STABILITY 

The Lyapunov exponents of a dynamical system characterize the rate of separation of pairs of initially adjacent 
trajectories [21;— 23; . Positive Lyapunov exponents provide a universal signature of Hamiltonian chaos. Their compu- 
tation requires a reliable algorithm for the integration of Hamilton's equations, which respects their symplectic nature. 
Unfortunately, integrating Hamilton's equations for the isochoric pentahedron has proven to be rather challenging, 
because the Hamiltonian is not analytically known in terms of the phase space variables. All derivatives of the Hamil- 
tonian with respect to the phase space variables must be computed numerically. This introduces numerical errors 
into the Hamiltonian equations of motion, which may not be relevant over short segments of the trajectory, but can 
introduce large errors in long trajectories, which are needed for the determination of the global Lyapunov exponents. 
To avoid this complication, we here pursue a simpler method which gives useful insights into the integrability and 
stability of the pentahedron as a dynamical system. 

The Kolmogorov-Sinai (KS) entropy gives the rate at which a dynamical system destroys information (creates 
entropy); it is given by the sum over all positive global Lyapunov exponents of the system. Lacking a reliable numerical 
method for calculating the GLE's, we will estimate the KS entropy in terms of the local Lyapunov exponents. In 
doing so, it is worth noting that local instability does not guarantee the global instability of a trajectory. For example, 
the well known x 2 y 2 system is almost everywhere locally unstable, but it admits nontrivial periodic trajectories that 
are globally stable. In general, the phase space averaged LLE's will provide an upper bound for the KS entropy, 
because the local divergence between two trajectories can be balanced by local convergence in other places, leading 
to a smaller global rate of divergence along the entire trajectory. On the other hand, while a small number of globally 
stable trajectories cannot be excluded in a dynamical system that is locally unstable in a large fraction of phase space, 
most trajectories will be globally unstable. 

To calculate the local Lyapunov exponents (LLE) of a Hamiltonian system we consider the small local deviation of 
a pair of trajectories 5z(t), which can be computed by linearizing Hamilton's equations: 



6z(t) = H{t,z)5z(t), 



(23) 
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where % is the Jacobian matrix of the Hamiltonian evaluated at z 



™'- :): r|# S) 



Solving (231 gives Sz(t) = exp(X z t)u where X z is an eigenvalue of %(z) and u is the associated eigenvector. Eigenvalues 
with positive (negative) real components are associated with exponentially diverging (converging) trajectories, purely 
imaginary eigenvalues represent a periodic motion. Although the LLE's of a system typically do not correspond to 
the true or global Lyapunov exponents (GLE's) they serve as a local measure of the departure of a pair of trajectories 
and give a good indication as to the stability of the system at a particular point. We define a point to be stable if 
all the positive real components of the eigenvalues of H are zero, i.e if max{5RA+} = 0. Our numerical explorations 
of the phase space have shown the Hamiltonian to be relatively smooth and continuous, computing the local stability 
over a finite grid should give a reasonable approximation to the overall stability of the system. 

If the entire phase space was found to be unstable it would be reasonable to conclude that the system is globally 
unstable and that the associated GLE's would be positive. However what can we conclude if we find a finite set of 
stable points amongst some set of unstable points? 

We compute the local stability over a grid of 2 20 points in the phase space, partial projections in the pi,P2 and 
qi,q2 planes showing the largest positive eigenvalues are shown in Fig. [18] and Fig.[l9j These figures show that most 
of the space is locally unstable (hot colors), however sizable stable regions (dark blue) do appear. The large values 
(light yellow) in the stability distributions shown are clearly correlated with the boundaries in the configuration space 



(see Fig. [16] and Fig. 17 1. Regions close to disallowed transitions appear very unstable. 

The distribution of stable regions is generally correlated with larger values of pi,P2 suggesting larger volumes. In 
Fig. [9] we show the distribution of volumes for stable and unstable regions, stable regions are much more likely to be 
associated with higher volumes. The mode of the stable curve is at V — 0.645, the mode of the unstable curve is 

r.l/4 

at V = 0.587. The volume of an equilateral triangular prism with unit face areas is = 0.658, stable regions are 
associated with more regular polyhedral configurations. As the figure shows this relationship is not cut and dried as 
the high mode of the unstable curve shows there is a sizable fraction of unstable configurations at relatively large 
volumes. 
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The volume distribution for almost stable regions < max{3?(A+)} < 1 is plotted against their stability in Fig 
There is a clear linear correlation between volume and stability in the highest density region of the figure. This 
correlation explains the relatively high peak of the inclusive and unstable curves in Fig.[9j these are associated with 
weakly unstable regions in the phase space. The majority of stable, or nearly stable, configurations correspond to 
those with a large volume. 

Knowing the local stability of a point is most useful if we have some measure of how likely it is that the a trajectory 
of the system will actually pass through this point. Ideally this measure would be computed by dividing the phase 
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FIG. 9: The distribution of volumes (black) as sampled over a grid of 2 20 points spanning the phase space. The volume 
spectrum of the stable regions is shown as the dashed red curve, the unstable regions are shown by the dashed green curve. 
Note that the area of all curves shown are normalized to one. 
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FIG. 10: The density of the positive real components of LLE's plotted against the volume of the system, the LLE's are restricted 
to the nearly stable region (0 < max{5ft(A+)} < 1 ). 



space into small volumes, integrating many trajectories, and then calculating the probability for a trajectory to reach 
each of these points. The Kolmogorov- Sinai entropy arises directly from this construction (cite). As this is currently 
numerically untenable we seek an alternative measure for the probability of the system being in a given configuration. 



A. Canonical Stability Estimates 

We consider the system as a member of a canonical ensemble, i.e we embed the system into a heat bath at some 
fixed temperature T = /3 . This non-isolated limit is a reasonable choice if our polyhedron is part of some extended 
interacting system of polyhedra. In the canonical ensemble we can assign the following weight to any particular 
configuration a of the system 

P(a) = UM-PH{o)) if H(a) > 
1 ' |0 if H{a) undefined, 1 ' 

where the second case explicitly assign zero weight to points where the Hamiltonian cannot be evaluated or where the 
numerical routines fail. The Kolmogorov-Sinai (KS) entropy for a system can be computed from the sum of positive 
GLE's, the largest positive LLE can serve as local estimate of this. The eigenvalues of H are opposite signed pairs, in 
unstable regions the secondary positive eigenvalue is usually small so this is a reasonable approximation. Using the 



canonical measure ( 25 1 we propose the following canonical instability measure as an estimate of the KS entropy 



%fT s _ Vxp(-/3ff(g))max{Sft(A+)}Ar 
[ >~ J u eM~l3H(a))da 1 ' 

where the region U is the set of unstable points in phase space, i. e. points where max{5R(A+)} > 0. In Table. [i] we 
compare estimates for the KS entropy along with the fractions of stable and unstable volumes of the phase space for 
the unit area triangular prism, the Henon-Heiles (HH) potential V(x,y) = \{x 2 + y 2 + 2x 2 y — 2/3y 3 ) and the two- 
dimensional classical Yang-Mills (CYM) potential V(x,y) — \x 2 y 2 . The latter are both well known non-integrable 
systems which display aspects of Hamiltonian chaos, see [3TJ [52 [3U] and references therein for more details. As the 
Table shows, the triangular prism has the largest relative volume of local stability, being almost twice that of the CYM 
and HH systems. The HH system is known to be non-integrable for larger energies, the CYM system is non-integrable 
for all energies. That our estimates of the KS entropy for the triangular prism is of a similar magnitude is strongly 
suggestive that it is also non-integrable and exhibits some degree of mixing. 



B. Microcanonical Stability Estimates 



To examine the energy (volume) dependence of the instability measure we consider the microcanonical ensemble. 
Estimates using the microcanonical ensemble are valid for isolated systems. If the system is ergodic then microcanon- 
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Name 


Momentum Range 


Position Range 


Stable Fraction 


Unstable Faction 


h(l) 


Henon-Heiles 


[-2,2] x [-2,2] 


[-2,2] x [-2,2] 


0.06 


0.94 


1.36 


Classical Yang Mills 


[-2,2] x [-2,2] 


[-2,2] x [-2,2] 


0.08 


0.921 


0.548 


Triangular Prism 


« [0,2] x [0,2] 


[0, 2tt] x [0, 2tt] 


0.13 


0.86 


0.38 



TABLE I: The fraction of stable and unstable regions along with the canonical stability measure ( 26 1 for the triangular 
prism system (with all faces set to 1), the Henon-Heiles potential and the two-dimensional classical Yang-Mills potential. The 
temperature is chosen as T — 1. 



ical averages are equivalent to averages over trajectories in the long time limit. The microcanonical density £l(E) 
is 



Q.(E) 



f 8{E - H{p,q)) dpdq, 
Jv 



(27) 



where the integral is computed over the entire phase space. The total volume distribution, the black curve, shown in 
Fig. [9] is proportional to the microcanonical density for the equal-areas triangular prism. 
The associated microcanonical proxy for the KS entropy is 



(hk(E)) 



n(E) 



max {5i(A+)} S(E - H(p, q)) dpdq. 



(28) 



In Fig. 11 we plot the microcanonical instability measure computed from LLE's for the unit area triangular prism and 
for a free particle in the Henon Heiles potential as function of the energy E. The instability measure of the triangular 
prism is high at small energies (pentahedron volumes) and decreases rapidly to zero at large energy. Small volume 
configurations exhibit a high degree of instability, because the system flops between the many available equivalent 
configurations. At larger volumes the set of possible volume preserving deformations is reduced and the instability 
decreases. As a comparison we have also computed the microcanonical instability measure for the HH system using 
(28). The HH system is stable for E < 0.1, and the instability rises rapidly for for larger energies. At low energies 



the system is known to be integrable 
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FIG. 11: Microcanonical instability measure hk(E), 
particle moving in the HH potential (b). 



computed using LLE's, for the equal area triangular prism (a) and a 



C. Intermediate Lyapunov Exponents 

Given the promising results of the local Lyapunov analysis we turn to examining intermediate Lyapunov exponents 
(ILE's). The LLE's suggest the existence of large areas of local instability however there may be some complicated 
dynamical conspiracy that allows for the existence of stable periodic trajectories. The ILE's provide a bridge between 
the LLE's, which are the zero time limit of the ILE's, and the global Lyapunov exponents (GLE's) which are the 
infinite time limit of the intermediate Lyapunov exponents. 

We compute the maximum positive ILE by the method of Benettin et al |31) . A reference trajectory is numerically 
integrated for some period of time z(t). A vector d is generated by some small deviation do from the starting point 
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of the reference trajectory which is then also integrated. The largest positive intermediate Lyapunov exponent after 
some time t n is 

A„ = -^— log^. (29) 
t„ — to do 

In practice the distances between the trajectories d(t) grow rapidly, to improve numerical stability the distances are 
periodically rescaled once they exceed some fixed threshold D. We define a rescaling 

«i = ||y {ti\d(t 1 )>D}, 



repeating this process as required the ILE ( 29 ) can be written as 



1 

A„ = V log a, , (30) 

t n ~ to %=1 

where <* = HI- 

The reference and deviation trajectories were integrated using an implicit symplectic Runge-Kutta (RK) integrator 
[32] . A computationally more expensive implicit method was chosen as the Hamiltonian cannot be readily separated 
into momentum and potential components making explicit RK splitting methods untenable. The implicit equations are 
solved iteratively. The iterations are initialized using an equistage method. Symplectic integrators explicitly preserve 
the two-form Q = dq A dp upon the phase space. A symplectic integrator generates maps which are themselves 
symplectic, orbits of the system under this map can be thought of as following a shadow solution to those of the 
exact Hamiltonian. Orbits under non-symplectic maps may follow the exact Hamiltonian more closely for small times 
but in general will diverge from the true dynamics over long time scales [321 133] . A fourth order Gauss collocation 
method was found to provide a satisfactory balance of accuracy and computational efficiency [32j [34] . The stability 
properties of these methods ensure that the dynamics does not depend upon the choice of time step for the integrator, 
the numerical convergence of the iterative process is slower at larger time steps. 

A large set of isochoric points were sampled from the Hamiltonian and used to generate trajectories and ILE's, this 
process was carried out for a set of volumes {0.35,0.45,0.55,0.65} so that the energy dependence of the ILE's could 
be computed. 



o V=0.35 □ V=0.55 
a V=0.45 x V=0.65 
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FIG. 12: The ensemble averaged intermediate Lyapunov exponents (ILE) computed for ensembles of trajectories with vol- 
ume {0.35,0.45,0.55,0.65}, the bars show standard errors. The dashed horizontal lines show the appropriate microcanonical 
estimates of the LLE's. 



The time evolution of the ensemble average of the ILE's for each of the volumes is shown in F ig.[L2| note that 
the y-axis is log scaled. The dashed lines show the microcanonical LLE values computed from (28 1. The ILE's 
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rapidly peak and then slowly saturate to stable long time values. Note that the rise time is different for the different 
volumes, the more unstable small volume (V = 0.35) ILE's grow faster than the more stable high volume case. In 
Fig. [13] the time evolution of the distribution of ILE's generated for the V = 0.45 case is shown. The distribution is 
centered around zero for small times, the mean converges to 0.2 by t — 40. The small time limit does not reproduce 
the microcanonical estimate of the LLE's. For very short times the LLE's have more information about the local 
instability of the system than the ILE's as they are generated from the full Jacobian of the Hamiltonian gradient. 
The LLE's tend to fluctuate around zero for small times as they are generated from random initial displacements 
which may be initially along a stable direction. As we observe in Fig. [12] over time the most unstable direction will 
dominate and the distributions appear to converge. 

The long time values should be reasonable approximations to the GLE's. The ILE's are larger at low volume and 
become very small at V = 0.65, this matches well with the microcanonical estimates. The unit pentahedron is strongly 
chaotic at small volumes and becomes more regular as the volume approaches it maximum limit. 



VI. CONCLUSIONS 



In our investigation of the phase space of the unit area triangular prism we have found a great deal of structure 
in the Hamiltonian and in the distribution of configurations. The phase space contains moderate regions of local 
stability and large regions of local dynamical instability. The distribution of local Lyapunov exponents appears to 
be correlated with the boundaries in the configuration space. We have calculated the average dynamical instability 
measures in the canonical and microcanonical ensembles and obtained values that are comparable to those found in 
well-known chaotic systems. The triangular prism differs from these and most other dynamical systems in that the 
stability increases with energy. Higher energy (high volume) triangular prisms are dynamically more stable than low 
energy (low volume) prisms. 

The large degree of dynamical instability found in our investigation of the isochoric pentahedron with unit area faces 
provides an encouraging starting point for a bottom-up investigation of the origin of thermal behavior of gravitational 
field configurations in loop quantum gravity. That the dynamical instability occurs in the simplest polyhedron where 
it can suggests that it will be a generic property of more complex polyhedra. Any coupling to other polyhedral 
configurations can be expected to enhance the degree of instability. We emphasize again the most interesting result 
of our investigation, that the rate of instability does not go to zero at small energy. In other words, the low energy 
dynamics of the isochoric pentahedron is not characterized by ballistic trajectories in phase space, but by rapid and 
energy independent diffusion among different degrees of freedom. At low energies, the pentahedron appears to be a 
fast scrambler of information. 

It is tempting to speculate whether this unusual property has implications for the quantum theory of gravity. 
Dynamical chaos in classical field theories is related to the instability of the perturbative effective action of the 
quantum field theory [351135]. Because the full effective action is known to be convex, and thus dynamically stable, 




FIG. 13: The distribution of intermediate Lyapunov exponents (ILE) for trajectories with V = 0.45 plotted at increasing times 
t — {0.1, 5, 10, 20, 40, 100}. The vertical red line shows the mean of the distribution. 
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for any quantum field theory [37], this instability is usually cured by some dynamical mechanism, such as mass 
generation through spontaneous symmetry breaking or confinement, or the breaking of translational invariance. None 
of these mechanisms are realized in the case of gravity, which is unique in being characterized by the simultaneous 
absence of an infrared mass scale and the existence of a ultraviolet scale, the Planck mass, which controls the quantum 
corrections to its classical limit. 

Returning to the dynamics of the pentahedron, we intend to report on the systematic numerical integration of a 
phase-space spanning set of trajectories in a future work. The preliminary efforts required to compute the ILE's have 
identified small regions of quasi-periodicity. The isochoric pentahedron thus does not exhibit full ergodicity. The 
isochoric pentahedron being a dynamical system with rich and fascinating dynamics, larger polyhedra with N > 5 
will likely be even more so. Investigation of the N — 6 system would be particularly interesting, as it would expose 
an additional aspect in the structure of the phase space as the hexahedron passes between the two dominant classes 
of six-faced polyhedra. 

The non-integrability of the classical triangular prism indicated by our analysis rules out a straightforward ap- 
plication of Bohr- Sommcrf eld quantization, which was successfully used to explore the quantum mechanics of the 
tetrahedron [19] . Methods from the study of quantum chaos such as the Gutzwiller trace formula |21) may provide 
some insight into the quantum behavior of the pentahedron and higher polyhedra. Variational methods or diagonal- 
ization in some suitable basis may also prove fruitful. 
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FIG. 14: Sections through the Hamiltonian in the pi,P2 space. The plots span a range of zero to n in qi, q^ with q^ increasing 
along the vertical axis and q\ along the horizontal axis. Cooler blues correspond to smaller volumes, brighter yellows correspond 
to larger volumes. 
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FIG. 15: Sections through the Hamiltonian in the qi, 52 space. The plots span the range [0, 2] in pi,P2 with P2 increasing along 
the vertical axis and pi along the horizontal axis. Cooler blues correspond to smaller volumes, brighter yellows correspond to 
larger volumes. 

FIG. 16: Slices through the phase space in the pi,P2 plane over gi,(?2 in [0,vr] showing the different configurations of the 
triangular prism. The projection here corresponds to that shown in Fig. |18| 
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FIG. 17: Slices through the phase space in the q\,qi plane over the full range of pi,P2, showing the different configurations of 
the triangular prism. The projection here corresponds to that shown in Fig.|19| 
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FIG. 18: The local Lyapunov exponents computed for sections in the pi,P2 over (71,52 in [0,7r]. The royal blue regions are 
stable and hot colors (red, orange, yellow) represent progressively more unstable regions. The code failed to converge in the 
light cyan regions. Note that here only LLE's in the range A £ [0, 1] are plotted here. 
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FIG. 19: The local Lyapunov exponents computed for sections in the ^1,92 over the full pi,P2 range. The royal blue regions 
are stable and hot colors (red, orange, yellow) represent progressively more unstable regions. The code failed to converge in 
the light cyan regions. Note that here only LLE's in the range A G [0, 1] are plotted here. 



